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Abstract 

This paper is devoted to existence and uniqueness results for classes 
of nonlinear diffusion equations (or systems) which may be viewed as 
regular perturbations of Wasserstein gradient flows. First, in the case 
where the drift is a gradient (in the physical space), we obtain exis¬ 
tence by a semi-implicit Jordan-Kinderlehrer-Otto scheme. Then, in 
the nonpotential case, we derive existence from a regularization proce¬ 
dure and parabolic energy estimates. We also address the uniqueness 
issue by a displacement convexity argument. 
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1 Introduction 

The continuity equation with a density-dependent drift 

dtp = div(pt>), with v = V[p] 

is ubiquitous in modeling and arises in a variety of domains such as biology, 
particle physics, population dynamics, crowd modelling, opinion formation... 
It should actually come as no surprise since it captures the dynamics of a 
population of particles following the ODE X = —v(t,X) where v = V[p\ 
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depends itself on the density in a way (local, nonlocal, attractive, repulsive 
etc..) depending on which phenomena (aggregation, diffusion...) one aims to 
capture and the type of applications. 

Of course, at this level of generality not much can be said on existence 
and uniqueness. There are however two cases which may be treated in a 
systematic way. The hrst one, is the regular case where V[p] is a smooth 
vector held whatever the probability measure p is, with some uniform bounds 
on some of its derivatives and p V[p] is Lipschitz in the Wasserstein 
metric. In this regular case, existence and uniqueness can be proved as a 
simple exercise by the method of characteristics and a suitable hxed point 
argument. This regular case (a typical example being that of a convolution) 
is however rather restrictive and for instance rules out diffusion. The second 
case where there is a general theory is the Wasserstein gradient flow case. 
In this case, at least at a formal level, v may be written as n = that 
is the gradient of the hrst variation of a functional S dehned on measures. 
In their seminal paper [10], Jordan, Kinderlehrer and Otto discovered that 
the heat how is the gradient how of the entropy functional 8{p) = f plog(p) 
which corresponds to the case v = ^. The theory of Wasserstein gradient 
hows has been very succesful in addressing a variety of nonlinear evolution 
equations such as the porous medium equation [IB], aggregation equations 
[3] or granular media equations [1]. This powerful theory is presented in 
a complete and detailed way in the reference book of Ambrosio, Gigli and 
Savare [1]. 

The purpose of the present paper is a contribution to the following general 
question: can one hope for an existence/uniqueness theory in the case where 
V is the sum of a Wasserstein gradient how term and a regular term (not 
necessarily a gradient). Our motivation for this question actually comes from 
systems. For instance, a simple but natural model, for the evolution of two 
(say) interacting species is: 

(dtpi = niApi + div(piV(F*pi + G*p2)), 

\dtP2 = J^2^p2 + div(p2V(iJ *pi + K-k P 2 )). 

When 1^1 = z /2 = 0 i.e. without dihusion, this is exactly the system studied 
by Di Francesco and Fagioli [8]. As emphasized in [8], if cross-interactions 
are symmetric i.e. G = K (or more generally G and K are proportional), 
this system has a (product) Wasserstein gradient how structure but this 
is certainly a restrictive and often unrealistic assumption in applications. 
This is why Di Francesco and Fagioli, still taking advantage of the similarity 
with Wasserstein gradient hows used a semi-implicit scheme a la Jordan- 
Kinderlehrer-Otto to obtain existence and uniqueness results. In [8], there 
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is no diffusion but clearly the structure of the system belongs to the mixed 
case where drifts can be decomposed as the sum of a Wasserstein gradient 
and a regular term. Of course, the semi-implicit scheme only makes sense 
when drifts are gradients. 

Regarding systems with a gradient structure and in the presence of nonlin¬ 
ear diffusion, our hrst contribution is to establish strong enough convergence 
at the level of the semi-implicit scheme to recover a solution of the PDE at 
the limit. The delicate step is of course to pass to the limit in the nonlinear 
diffusion term, which can be done thanks to the powerful flow interchange 
argument of Matthes, McCann and Savare [13] in a similar way as in the 
work of Di Francesco and Matthes [9] . We will then address the nonpoten¬ 
tial case in which the drift may contain a nongradient (but regular) part. 
This case cannot be attacked by the semi-implicit minimization scheme and 
we will prove existence by suitably regularizing the diffusion and using stan¬ 
dard parabolic energy estimates. Finally, we will derive an uniqueness result 
from displacement convexity of the energy. 

The paper is organized as follows. In section O we consider the potential 
case, introduce a semi-implicit scheme a la Jordan-Kinderlehrer-Otto [101 
and state a hrst existence result. Section [3] is devoted to the proof of this 
existence result. Section 0] extends the result to the case of systems (again 
in the case where all drifts are gradients). Section [5] proves existence for 
the non-potential case. The hnal section [6] shows uniqueness by a simple 
displacement convexity argument. 

2 The potential case and the semi-implicit 
JKO scheme 

Our aim is to solve the following nonlinear diffusion equation: 

dtp = div(^pV(F'(p) + f7[p])), p\t=o = Po, (2.1) 

on (0, -|-cxd) X where T'^ := denotes the hat torus (we take periodic 

boundary conditions to simplify the exposition, we refer to the work of the 
second author for extensions to or a bounded domain) which we 
identify with the unit cube [0,1]'^ equipped with the quotient distance: 

y) := inf^ \x — y + k\^. 

Denoting by P(T‘^) the set of Borel probability measures on T'^, we assume 
the following assumption on the map p G V{T^) i-G- U{p)\ 

'ip G P(T'^), U[p\ G and U[p] > 0, (2.2) 
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(2.3) 


sup {||Vf/[p]||Loo + ||(A[/[p])+||loo} < +CX) 


pGP(Td) 


and there exists a constant C such that for all (p, v) G P(T'^) x ViT'^) 


||Vf/[p] - V17 HIU.o(t.) < Chh2(p,^), 


(2.4) 


with W 2 {p-, v) denoting the 2-Wasserstein distance between p and v i.e. 



where n(/i, z/) is the set of transport plans between p and u i.e. the set of 
Borel probability measures on T'^ x having p and u as marginals. It is well 
known, see [231 El Eo], that W 2 metrizes the weak star topology on V(T‘^) 
so that {V(T^), W 2 ) is a compact metric space. 

As for the nonlinear diffusion term div(pV(A''(p))) it is convenient to 
rewrite it as: 


div(pV(h;'(p))) = AF\p) 


where F'(t) := tE'{t) — E(t) so that F"(t) = The typical energies E 

we have in mind are the following classical examples 

• Entropy: E{t) := flog(t) so that F'{t) = f, F"{t) = 1 (which thus gives 
a linear diffusion driven by the laplacian), 

• Porous media E{t) = with m > 1 so that F'{t) = {m — 1)P^, 

m = 

We shall assume that E is a continuous convex function on M_|_ with 
E{0) = 0, E is of class on (0, + 00 ) and that there are constants C > 0 
and m > 1 such that 


E"{t) > — F'{t) = tE'{t) - E{t) < C{1 + r), Vf e (0, +CX)). (2.5) 

(_y 


Of course, these assumptions are satished in the examples above correspond¬ 
ing respectively to linear diffusion and the porous medium equation. 

Finally, we assume that the initial condition po ^ ViT^^) satishes 


/ E{pq{x)) dx < -fcxo 


( 2 . 6 ) 


which with fl2.5p in particular implies that po G L™'(T°') and E'(po) G (¥“*). 
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A weak solution of fl2.ip then is a curve t G (0, +oo) i-)- p{t,.) G P(T‘^) 
such that F'{p) G L[q^(( 0 , +cxo) x T'^) and 


'•+00 


Ijd 


{dt4>p + A(f>F'{p) — VU[p] ■ V(j)p) dx dt = — j 0(0, x)pq{x) da; 

(2.7) 


/Td 


for every 0 G Cc ([0; + 00 ) x T^). 


Theorem 2.1. Assume 
least one weak solution. 


■H2.A)-^2.4^-^2.5\)-H2.6]). then H2.1\) admits at 


The complete proof of this result is given in section [3l This proof is 
strongly based on a semi-implicit version of the Jordan-Kinderlehrer-Otto 
(JKO) scheme [10] as in Di Francesco and Fagioli [S]. More precisely given a 
time step h > 0, we construct inductively a sequence p^ G P(T^) by setting 

= Po and, given we select as a solution of 

Afr-I + «('>) (2-8) 

where 

S{p) := IF := f UMdp. (2.9) 

1-|-cx), otherwise, Jjd 

Note that assumption fl2.5l] ensures that S controls from above p"^dx if 
m > 1 and plog(p)dx if m = 1, so in any case sublevels of are weakly 
relatively compact in L™'(T'^). 

By standard lower semicontinuity and compactness arguments, it is clear 
that fl2.8p possesses solutions so that one can indeed generate a sequence 
(Ph)kGN by the semi-implicit JKO scheme 02.81) . It is even uniquely dehned 
(but we won’t really need it in the sequel) because each p^ remains absolutely 
continuous with respect to Lebesgue’s measure, so that p 1 —)■ 110^(p,p^) is in 
fact strictly convex (see Proposition 7.19 of [20|) and the other terms £ 
and 7/(.|p^) are convex. We hnally extend in a piecewise constant way the 
sequence (p^)fcgH i.e. set: 

Ph(t ,.) := p\ for t G {{k — l)h, kh] and /c G K. (2-10) 

The proof detailed in the next section consists in showing that as h —?■ 0, one 
may recover a limit p which satishes 02. ip . This is the same strategy as in [8] 
but the tricky part consists in passing to the limit in the nonlinear diffusion 
term F'{ph). This will be done thanks to the powerful flow interchange 
argument of Matthes, McCann and Savare [T3| in a similar way as in the 
work of Di Francesco and Matthes [9] . 
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3 Existence proof 

The proof of Theorem 12.11 is divided into three parts. The hrst two parts 
concern a priori estimates on ph and the last part consists in showing that 
passing to the limit in the Euler-Lagrange eqnation of fl2.8p actnally enables 
us to recover a solution of fl2.ll) . The discussion on uniqueness is deferred 
to the hnal section O Of course, it is enough to work on a hxed hnite 
time interval (0,T), which we shall implicitly do below, we thus also set 
N := \^] + 1. In what follows C (respectively Ct) is a generic (resp. only 
depending on T) constant whose value may vary from one line to another 

3.1 Basic a priori estimates 

From the very dehnition of the JKO semi-implicit scheme 02.81) we have for 
every k\ 

Tn'-|(p;;+‘, pI) < £(pt) - f (p^) + £ U{pi]dU - pI*')- (3.1) 

Recall then that the 1-Wasserstein distance Wi is dehned by: 

Wi{p,u)-.= inf { / d(x,|/)d 7 (x, 2 /)|, 

7en(p,i/) i. jfdy^fd ) 

so that by Cauchy-Schwarz’s inequality 


Wy{p,v)<W2{p,v). 

The well-known Kantorovich duality (see [2311211EQ]) states that W\ can be 
also be expressed as 


lKi(p, v) = sup < / 0d(p — u), (j) 1-Lipschitz on 

so that 0d(p — z/) is less than the Lipschitz constant of 0 times lTi(p, zz). 
Thanks to these considerations, assumption fl2.3p and Young’s inequality, we 
get 

j^^( 7 |p|]d(p‘-p!;+') <civ2(p!;+',p‘) < Tn^|(pyi,p|) + c2ft. (3.2) 

Together with 03.ip . this gives 


inw. p‘) < f(rf) - £(pr)+c/i 


4ft 




(3.3) 
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summing between 0 and N and using the fact that E is bounded from below 
and fl2.6p gives 


jr E Pk) S f(Po) - e(Pk) +CNh<C{l+T), (3.4) 

^ k=0 

as well as 

£(pi)<£(po)+Ckh, (3.5) 

which, thanks to fl2.6p and fl2.5p . also gives 

sup \\phit,.)\\Lm<C{l + T)iim>l, (3.6) 

telo,T] 

and 

sup [ ph{t,x)\og{ph{t,x))dx < C{1+ T), a m = 1. (3.7) 

tG[0,T] JT'* 

With fl3.4p . we also have the Holder like estimate: 

W 2 {ph{t, .),Pk{s, .)) < CV(1 + T)^\t -s\ + h, V(s,f) e [0,T]2. (3.8) 


Using fl3.8p and rehned versions of Ascoli-Arzela Theorem (see [T]) and fl3.6p - 
fl3.7lh one deduces the existence of a vanishing sequence —)■ 0 and of a 

p e (P(T'^),W 2 )) nL°°((0,T),L”^(T'^)) such that 

P/in ^ P in 7.”"((0,T) X T'^),!™ sup lU 2 (p/x„(f,.), p(f,.)) = 0. (3.9) 

" t£[0,T] 

Now, using fl2.4p . we deduce that 'VU[ph^] converges to VU[p] in L°°((0, T) x 
T'^) and then 

PhyU[p^J - pVU[p] in L™((0,T) X T"). (3.10) 


3.2 Refined a priori estimates by flow interchange 

This is the key step in the proof which enables us to obtain strong con¬ 
vergence, in what follows we essentially follow similar arguments as in Di 
Francesco and Matthes [9]. For u G P(T‘^), set 


Ent(i/) := 


z/((r) log(z/((r))dx, if iy\og{iy) G 
-|-oo, otherwise. 


For u G V{T‘^) with Ent(i/) < -|-cx) let us denote by := the 

solution at time t of the heat equation: 

dtp = Ap, p |4=o= (3.11) 
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It is well-known since the seminal work of n that the heat flow can be 
viewed as the gradient flow of Ent for W 2 (see [T] for the theory of gradient 
flows in metric spaces). Moreover the fact that Ent is displacement con?;e20, 
gives (see [1] Theorem 11.1.4, [6], [IS]), the following evolntion variational 
eqnality: 


1 ci+, 

2 dt 


p) < Ent(/i) - Ent(e‘^i^), Vf > 0, p G V{T) (3.12) 




where we have used the notation: 


at s—> 0 + ® 


Taking as a competitor in the minimization fl2.8p gives 


0 < 


1 

2h dt 


Wi{e-pi 


tA k+1 



(3.13) 


Since is smooth for f > 0, we can directly compute: 



which, with fl2.5p gives that for some positive constant A > 0 


We then have 

- ^(f > Alimmf j(‘ |V((e‘*^py‘)t)|Mxd(. (3.14) 

^See section [5] for a precise definition. Here, we are working on T”*, but we should not 
worry about it, it is just if we were working on with periodic functions only. The optimal 
transport map between absolutely continuous periodic measures is well-known, it is given 

I 1^ 

by the gradient of a convex function F such that F{x) — is periodic (see Cordero- 
Erausquin 0) and which is characterized by a Monge-Ampere equation. Displacement 
convexity of the entropy on the can therefore be proved as in the euclidean case. 
Another way to see this is to remark that Bochner’s formula on is just the same as 
in this does not change the Ricci curvature and thus, thanks to a celebrated result 
of Lott and Villani m and Sturm (see in particular the proof of Theorem 4.9), this 
does not change the displacement convexity of Ent (with respect to Lebesgue’s measure). 






In a similar way, for f > 0, we have 


^(M(e“p!;+'|p|)) = A(l7|p‘])e“p!;+‘ 

and the right hand side is nniformly bounded from above thanks to 02.31) . 
With 03.12^ - 03.13^ - 03.141) this gives: 

Ahliminf [ [ )|Mxdt < (Ent(p^) - Ent(p^+^)) + Ch. 

Since converges strongly to in L™' as s ^ 0+, (e^^p^'*'^)^ con¬ 
verges strongly to in By lower semicontinuity we deduce that 

V(p^+^)f e L2 and 

h[ |V(py')¥)|Mi<C(Ent(pE)-Ent(py') + fc). (3.15) 

JT'i 

Summing from /c = 0 to iV — 1 and using the fact that Ent(po) is hnite gives 

I [ \V{ph)'^)fdxdt<CNh + CiEnt{po)-Ent{p^))<C{l + T) 

Jo 

which, with 03.61) . also gives 

||p/i(t, .) 2 — ^T- (3.16) 

We then observe that since the injection of in L^(T'^) is compact 

and since p i—>■ maps continuously L‘^(T'^) into L'^(T'^), sublevel sets 

of p 1-4 \\ph(t, .)~\\Hi{Td) are strongly relatively compact in L™(T^). Now 
arguing as in Di Francesco and Matthes [9] i.e. using the rehned version of 
Aubin-Lions Lemma provided by Theorem 2 of Rossi and Savare [TH] gives 
that the family {ph)h is relatively compact in L™((0, T) xT^). The conclusion 
of this step is that 03.9p can be strenghtened to 

Ph„ P strongly in L™((0,T) x T'^). (3-17) 

Now, thanks to the second part of 02.51) and Krasnoselskii’s Theorem (see [7], 
chapter 2) p i-4 F'{p) is continuous from to and then 03.17p implies 
that 

F'iphrJ) F\p) strongly in L^((0,T) x T'^). (3.18) 
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3.3 Passing to the limit in the Euler-Lagrange equa¬ 
tion 


Now, we write the Euler-Lagrange equation for 02.81) as in the seminal work 
of Jordan, Kinderlehrer and Otto ra for the Fokker-Planck equation and 
Otto im for nonlinear diffusions. Let be a C°° vector held on and 
denote by Xt the how of 

= ^(Xt(a:)), Xq{x) = x. 

Dehne then := Xt^p^^ so that the change of variables formula gives 

pl+^ = Vt{X^)dei{DXt). (3.19) 

Since p^^ solves 02 . 8 |) . we have 


1 

0 <- 

2h dt 


^ / \t=0 




(3.20) 


Let then 7 ^ G n(p^,p^''“^) be such that 




Choosing k{x, y) G such that d^i^x, y) = \y + k{x, y) — xp, by dehnition of 
Vti we then have 


Wi{uupi)< [ \Xt{y) + k{x,y) - x\‘^d-f^{x,y) 


and since Xt{y) = y + t^{y) + o{t), we get 


1 (i+ 
2h dt 


^ / \t=0 


fjdy^jd 


Uv) ■ (^^-j^)d7h(x,P)- 


(3.21) 


As for the diherentiation of S{vt), following [in]-[IZ], using 03.19P we write 

observing that for p > 0 and a > 0 
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then thanks to fl2.5p . fl3.6p . the fact that det(Xt) = 1 +Miv(,^) + o{t) (with a 
uniform o{t)) and Lebesgue’s dominated convergence Theorem, one obtains: 

d+ 


In a similar way, 


dt 




= - / 


|t=0 


I'Id 


k-\-l 

'h ■ 


(3.22) 


(3.23) 


Combining fl3.2ip - fl3.22p - fl3.23p and fl3.20p . and applying the previous to both 
^ and —^ gives the Euler-Lagrange equation 


/T'^xT'^ 


i(y)-(x-y)A',^{x,y) = -h F'(pr)div«) + ft / V(7|p‘] ■ f 

(3.24) 

for every smooth vector held Now let 0 G C'^([0,T) x T'^) (which we 
extend by 0(0,.) on (—h, 0)), we then have 

N 

dt^Ph = X] / .) - 0((/c - l)h, .))pI 


f 

/o 


k=0 

N 


/'Id 


^2 [ Hkh,.){Ph-p’t^) - [ 0(0, .)po 


k=0 

N 

k=0 ■ 


'T^xT^ 


{(l){kh,x) - (j){kh,y))d-f^{x,y) - / 0(0, .)po- 


ITd- 


Using a second order Taylor-Lagrange formula gives 


J'fd^'fd 

with 


{(j){kh,x)-(l){kh,y))d-f^{x,y) = 


I'fdy^'^d 

IfiftI < \\D^/>\\L~WM+\pi). 


V0(/ch, y) ■ {x-y)d-f^{x, y) +RI 


With fl3.3p this gives that X]fc=o \^h\ — Crh, so that applying fl3.24p to 
^ = V(j){kh ,.), and using the fact that, with fl2.4p . Cauchy-Schwarz inequality 
and ([3ZD, Ef=o II“ Vf/[p^+^]||Loc < Ct, we hnally get 


'0 JTd, 


dtfpph = 


>0 JTd 

f 


F\ph)^4> 


\/U[ph] ■ V(j)ph - / 0(0, .)po + Sh 

Jfd 


where 6h goes to zero as h —?■ 0. Thanks to fl3.10p . fl3.17p and fl3.18l) we may 
pass to the limit on the vanishing sequence to hnd that the limit p is a 
solution of fl2.1lh This completes the proof of Theorem 12.11 
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4 Extension to systems 

Let us now consider the extension of fl2.1j) to systems for the evolution of I 
densities p := (pi,..., p;) of interacting species: 


dtpi = + Ui[p\)) = 0, pi|t=o = A,o, f = 1, • • •,/, (4.1) 

on (0, +oo) X T'^. Let us assume that for every i = 1,..., / the map p G 
P(T^)' ^ Ui{p) fullhlls 

Vp G V{T^)\ Ui[p\ G and Ui[p\ > 0, (4.2) 

sup {||Vf/i[p]||i,-+ ||(Af/i[p])+||i,oo} <+CX) (4.3) 

and there exists a constant C such that for all (p, v) = ((pi,..., p;), (i^i,..., vi)) G 
P(T'^)^ X P(T'^)^ and every i 

i 

\\VUM - Vf/,H|Uoo(T.) <Cj2W2{Pj,iyj). (4.4) 

1=1 


As in the previous section, we assume that for each i, is a continuous 
convex function on M_|_ with L^i(O) = 0, Ei is of class on (0, +cxo) and that 
there are constants C > 0 and > 1 such that 

fmi-2 

E':{t) > tE[{t) - E,{t) < C(1 + r^), Vt G (0, +oo). (4.5) 

Finally we assume that the initial condition po G P(T'^) satishes 

[ Ei{pi^o{x)) dx < +00 (4.6) 

*=i 

which with fl4.5p in particular implies that pi^ G L”^*(T'^). The semi-implicit 
JKO scheme then takes the following form: given a time step h > 0, we 
construct inductively a sequence p^ G ViT'^Y by setting p° = po and, given 
p^ G V{T‘^y we select p^^ as a solution of 

t 4,.) + +«(PI4)} (4.7) 

2=1 

where 

ELi Ijd Ei{pi{x)) dx, if Ei{pi) G L^, 

I -|-oo, otherwise. 
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and 


i=l 

Extending in a piecewise constant way the sequence {p’^)km defines the 
P(T'^)hvalued curve: 

Ph{t, .) := Pk for t G ((/c — l)h, kh] and /c G N. (4.8) 


Arguing exactly as in the proof detailed in section [21 there is strong conver¬ 
gence in ((0, T) X T'^) of a sequence ph„ to some limit curve p and 

passing to the limit in the Euler-Lagrange for fl4.7p exactly gives: 


Theorem 4.1. Assume then ( j^.ij ) admits at 

least one weak solution. 


5 The non potential case 

We are now interested in the case where the drift may not be a gradient. 
More precisely, we consider the following nonlinear diffusion equation: 

dtp - div(pV(E'(p))) -h div{pV[p]) = 0, p\t=o = po, (5.1) 

on (0,T) X T'^. Denoting by the dual of we assume the 

following regularity on the drift term V[p\. 

Vp G r\V{T^),V[p] G sup {\\V[p]\\l^ + ||VE[p]|Uoo} < +oo 

peL^nV{Td) 

(5.2) 

and for every i? > 0, there exists a modulus uji such that, for every (p, u) G 
(L^(T“') n'P(T'^))^ such that ||p||j|/-i('j'd) < R and < R, one has 

||l/[p] — E[z/]||j;^2(Td) < l^fidlp — ^^||_ff-l(T‘i))- (5-3) 

Typical examples of velocity fields p i-G- V[p\ that satisfy the above assump¬ 
tions are those of the form V[p\{x) = B{x,y)p{y)dy with B smooth 
enough (but not necessarily a gradient with respect to x). 

As before, E is convex on R_|_ and we define F'{t) := tE'{t) — E{t) so 
that E"(t) = tE"{t). We make the following assumptions on E (which are 
satisfied for instance when E{t) = t™' with m > 1 or E{t) = tlog(t)): 

E G C^{R+, M), F(0) = F'(0) = 0, F is convex, (5.4) 

E” is nondecreasing, and for every p > 0, F"(p) > 0 (5-5) 
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and there is a constant C > 0 such that 


F’{p) < C(1 + + F(p)), Vp e R+. (5.6) 

As for the initial condition po we assume that it is a probability density 
such that 

Po e H\Y), Fipo) e Li(T"), F'(po) e H\r). (5.7) 

A nonnegative weak solution of the PDE 

dtp - A(F'(p)) + div(pl/[p]) = 0, p|t=o = Po- (5.8) 

is by dehnition a function p G L^((0,T) x T'^,M_|_) such that 

F'(p) G L\{0,T),H\T‘^)) 


and 

I {—dt(j)p + VF'(p) • V(/) — pld[p] ■ V(/)) dxdt = / 0(0, x)po(x)da; 

-j-d j-fd 

(5.9) 

for every 0 G C^([0,T] x T'^) such that 0(T,.) = 0. 

Before we proceed to the existence proof, we need some preliminary re¬ 
sults. Let us hrst study the continuity of the drift term p = p(t, x) i-G- 
V[p{t, •)](a^)- If is easy to see that when 05.21) and 05.3p are satished and p" 
converges strongly in L^((0,T) x T'^) (hence in L^((0, T), i7”^(T'^)) to some 
p then V[p^] converges to V[p] in L^((0,T) x T'^), but we wil need a variant 
in the sequel: 

Lemma 5.1. Assume that 115. gl) and hi. A) are satisfied. Let p” be a sequence 
in L^((0,T) X T'^) such that dtp^ G L^((0, T),//“^(T'^)) with 



sup ||(9ip"||^2((o,T),ir-i(T‘i)) <+ 00 , (5.10) 

n 

and p G L^((0,T) x T'^) such that p'^ ^ p in L^((0,T) x T'^), then V[p'^] 
converges to V[p] strongly in L^((0,T) x T^). 

Proof. First observe that O5.10p implies that for some constant C one has 
||p"(t,.) -p”(s, .)||^^-i < C^/\t - s\, Vn, V(t,s) G (0,T)^ (5.11) 


Let t G (0,T) and for h G (0,f) dehne 


Pt,h{x) ■= 


h 


't—h 


p^{s,x)ds, Pt^y^ : = 


h 


p(s, x)ds 


' t—h 
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thanks to 05.111) . we obtain, for every n, t and h\ 


mt, ■) - < cTa, ut ,.) - < cVTl (5.12) 

For fixed h > 0, in L^(T'^) as n —)■ cxo, and since the imbedding 

of L^(T‘^) into is compact we also have W'f^h ~ 'Pt,h\\ 

n ^ 00 . We then get 

Wp^i^i ■) — P{^i Ollr^-i < 2Cy/h + Wp'l,^ - Pt,h\\H-^{T:d) 

from which we deduce that ||p"'(t,.) — pit, OIIh-i tends to 0. Thanks to 05.31) . 
this implies that ||l/[p"'(t,.)] — V[pit, .)]\\i 2 (jd) tends to 0. The claimed 
convergence then follows from 05.2p and Lebesgue’s dominated convergence 
Theorem. 

□ 

We now introduce a regularized nonlinearity to approximate 05.Sp by a 
uniformly parabolic equation as follows. Let e G (0,1), let 6^ and be 
respectively the smallest p for which F"ip) > e and the largest p for which 
F''ip) < e~^. Let then F^ be dehned by 

r F(5,) + F'i5,)ip - 5,) + I(p - 4)2 if p e [0,4]; 

F4p):=<^F(p)ifpe[4,M,], (5.13) 

[f(M 4 + F'(M4(p - M,) + i(p - M,f if p > M,. 

Clearly, by construction is convex and on M+ with 

e < F'J < ^ on M+ (5.14) 

and Fg converges pointwise to F since 4 and converge respectively to 0 
and +CXD. In fact, this approximation also has good T-convergence properties: 

Lemma 5.2. Let 9 G L^((0,T) x T'^,M_|_), then 

I F{6{t,x))dxdt (5.15) 

-jd 

moreover if 9^ G L^((0,T) x T‘^),M+) weakly converges to 9 in E L^((0,T) x 
T'^), then 


lim 

£->0+ 


Ijd 


Fei9it, x))dxdt = 


liminf f f F^i9eit, x))dxdt > 

Jo 




F(6'(t, x))dxdt 


(5.16) 
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Proof. Fatou’s lemma first yields 

lim inf / / 

£->■ 0 + Jq Jjd 

on the other hand 


F^{9(t,x))dxdt > / / F{9(t,x))dxdt 


'0 JT'i 


'0 Jt<^ 


F^{9{t,x))dxdt < / / F{9{t,x))dxdt+ 


'0 Jfd 



{F,{9)-F{9))dxdt 




since the second term in the right hand side converges to 0 , we easily deduce 
fl5.15p . Let us now assume that 9e € L^((0, T) x M+) weakly converges to 
6 * in G L^((0,T) x T*^). Let 7 > 0 (fixed for the moment) and denote by F'^ 
the function dehned by 


F\p) 


F{p) if p e [ 0 , 7 ], 

F( 7 ) + F'( 7 )(p- 7 ) if p > 7 


by construction F"^ is convex and below F. For e > 0 small enough so that 
7 G [(fejMg], we similarly define 


F]{p) 


Fe{p) if P G [ 0 , 7 ], 

F(7) + F'(7)(p-7) if p > 7 


so that Ff is convex and coincides with F'^ on + 00 ). We then have 


lim inf / / 

e-i> 0 + Jq Jjd 

[T - 

> lim inf / 


Fe{9e{t.,x))dxdt > lim inf 


£—)- 0 ^" J Q 


/ F'^(9i;(t,x))dxdt+ \imini 
e^0+ Jq Jjd £—>-0 + 



Ff{9^{t, x))da:df 

{FM) - m)) 


{6e<&e} 


the second term converges to 0 whereas by weak lower semi-continuity (thanks 
to the convexity of F'^) we have 


hence 


liminf f f F^{9s{t,x))dxdt > 
£->■ 0 + Jq Jjd 




F'^{9(t, a;))da;dt, 


lim inf 

£^ 0 + 




Fe{9s{t,x))dxdt > sup 

7>0 



F'^{9(t, x))dxdt 


and then fIS.lbp easily follows from the previous inequality, the fact that 
F'’' converges monotonically to F and Beppo-Levi’s monotone convergence 
Theorem. □ 
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Theorem 5.3. Assume H5. 7| j, then / (5.g|) admits 
at least one weak nonnegative solution. 


Proof. The proof proceeds in three steps. 

Step 1: Regularized equation. We hrst prove existence of a weak solution 
to the regularized equation: 


dtp^ - A{F'{p^)) + div(pW[p^]) = 0, p%=o = Po- (5.17) 

For hxed u G L^((0,T) x the linear parabolic equation; 

dtp’' - div(F"(i^)Vp^) + div(pW[i^]) = 0, p'\t=o = po (5.18) 

admits a unique weak solution which we denote := moreover p^ 

is nonnegative by the maximum principle and G L^((0, T), i7^(T'^)) fl 
C([0, T], L^(T'^)) and dtp G L^((0, T), i7”^(T'^)) and more precisely, thanks 
to fl5.2p and fl5.14p there is a constant such that: 


[ [ (iVp*"!^ + (p'')^)da;dt + /" \\dtp'\\H-i < Ce. (5.19) 
Jo Jt^, Jq 

Thanks to 05.31) . it is easy to check that is a continuous map of T^((0, T) x 
T'^,M+). In addition, 05.191) and the Aubin-Lions lemma (see [2], [2T]) im¬ 
ply that T^(L^((0,T) x T'^,M+)) is relatively compact in L^((0,T) x T'^). 
Schauder’s hxed point Theorem then ensures that admits at least one 
hxed point p^ i.e. a solution of 05.171) . 

Step 2 : A priori estimates. We aim now to derive estimates independent 
of £ on p^. Let <5 > 0 such that 5 G {6s, M^), we then take (p^ — 5)+ as 
test-function in 05.17p (which is actually licit since this test-function belongs 
to L^((0, T), i7^(T'^))) integrating between 0 and t G [0,T] this yields 


{p’'—6)+),H^ds+ 



0 J{p^>S} 


F''(/)|Vpl2 = 



pW[p"]-Vp" 


0 J{p^>S} 


hence, using Young’s inequality 


^\mt,.)-6u\\i.--\\{po-6u\\i.+ 



F''{p')\Vp'\^ 




< C 



|Vp^P + 


0 J {p^>S} 
1 


0 J{p^>S} 





|VpV 

lo J{p^>S} ^ 


2u 


0 j{p^>&} 



(p- 


e\2 


0 J{p’^>S} 



!(/ - + ^^l) 
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since F"{6) = F"{6) > 0 and F" nondecreasing, we can choose z/ small 
enough so that the hrst term in the right hand side is absorbed by the left 
hand side of the inequality. Gronwall’s lemma then gives 

sup ||p^(t, .)||l 2 < G (5.20) 

tG(0,T) 

for a constant C that does not depend on e. Next we take F^{p^) as test- 
function which similarly gives: 


Ijd 





F,(p‘(t,.))- I F,(p„)+ I I \VF:(p’W= ( / p‘V[p‘]-VF:,(p‘) 


0 JTd, 


<G(’^ 



0 

|vq(/)P + ^ 

Iq Jjd zu 




ip^ 


0 Jt'^ 


using 05.201) and chosing u small enough we thus get 


sup / F,{pe{t,.)) + 

ZG[0,T] Jt^ 


'0 Jt‘‘ 


vf;(p‘)i^ < c 


(5.21) 


for a constant C not depending on £. Next we use 05.61) and 05.20p - 05.2ip to 
deduce that 

sup [ Fl{p^) < C (5.22) 

te[o,r] Jjd 

together with PoincarGWirtinger inequality, using again 05.2ip . this gives 


(5.23) 


(5.24) 


l|-^e(p‘^)llL2((0,r),Hi(T'i)) < G. 

Step 3: Passing to the limit. Let us set 

n" ;= F'{p^), a" :=Vu^ - p^V[p^] 
so that 05.171) can be rewritten as 

dtp^ = - div(p^G[p^]) = div(a^), p%=, = po. (5.25) 

We know from the previous step that 

IIP^I|l°°(( 0,T),L2(T'Z)) + ||o'^||L2((0,T),L2(Td)) + \\u^\\L^{{0,T),m(Td)) < C (5.26) 
as well as 


ll^4P^IU2((0,T),ir-l(T‘Z)) < G. 


(5.27) 
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Passing to subsequences if necessary, we may therefore assume that 
p" ^ p in L2 ((o, T) X T^*), ^ m in T), H\T^)) (5.28) 

and thanks to Lemma 15.11 and 05.271) , we have 

^ a :=Vu- pV[p] in ^^((o, T) x T^). (5.29) 

Obviously one then has: 

dtp = div(a) = Au- div(pl/[p]), p\t=o = po. (5.30) 

So to establish that p is a weak nonnegative solution of 05.8p . it is enough to 
prove that u = F'[p). Thanks to the convexity of F this amounts to prove 
that 


'o 


F{9{t, x))dxdt > 


'0 JTd, 


F{p{t, x))dxdt + 


'O JTd 


u{6 — p)dxdt 


for every 9 G L^((0,T) x T'^,M+). By dehnition of we know that 


(5.31) 


/o 


F^{9(t, x))dxdt > 


/o 


Feip^it, x))dxdt + 


/o Jt‘^ 


u^{9 — p^)da;dt. 

(5.32) 


Let us prove that 



For that purpose, let be the potential dehned by 


(5.33) 


-AV^" = p^ / ^|J^ = 0, e H\r). 


(5.34) 


Thanks to 05.201) . we have -0^ G L°°((0, T), iL^(T'^)) with a bound indepen- 
dendent of e: 

II V'0^||ioo((o,T),iyi(T‘^)) < C. (5.35) 

As for the time derivative of we observe that 

-A{dtip^) = dtp"" = div(cr'') 

so that, thanks to 05.261) . we have dt'V'i/j^ G L^((0, T) x T'^) and more precisely 
||5tV'0^||i2(('o_T)xT‘^) ^ I|<^^I|l2((0,T)xT‘*) ^ C 
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this proves that is bounded in iJ^((0,T) x T'^), hence converges in 
L^((0,T) X T'^), up to an extraction if necessary, to ^|J given by 


- AV^ = p, [ ^P = 0, 


(5.36) 


Weak convergence of and strong convergence of in then give 


lim 

£ 



fjd 


= lim 


^ Jo 




'o Jt'^ 


VmVV’ = 


/o Jt'* 


up 


which establishes fl5.33p . Next, we use Lemma [521 letting e tend to O'*", using 
fl5.33p we obtain inequality 05.3111 which proves that u = F'{p) and so p is a 
weak solution of 05.811 . concluding the proof. 


□ 


The previous arguments again clearly adapt to systems. More precisely, 
let us consider the system for the evolution of I densities p := (pi,..., p;): 

dtpi - A{F-{pi)) + div{piVi[p]) = 0, pi\t=o = Pi,o (5.37) 

on (0,+oo) X T'^. Assuming that each function A) satishes 05.411 - 05.5^ - 05.6p . 
that the initial conditions are probability densities which satisfy 

p,,0 e i7'(T"), F,(p,,o) e L\T^), T;'(p,,o) e i7'(T"), W = 1,...,/, (5.38) 


and for every f = 1,..., /, the maps Vi satisfy 

Vp e L\i:'^)\Vi[p] e w'’“(T^), sup {\\V[p]\\l^ + \\^y^[p]\\L^} < +oo 

peL2(Trf)* 

(5.39) 

and for every R> 0, there exists a modulus ur such that, for every (p, u) G 
^ 2 ('jd)Z ^ and ||z 2 ||j:^-i(Td)i < R, one has 


L 

— hiH||L2(T'i) < \\Pj ~ ^jWn-^iTd) 


(5.40) 


i=i 


A direct adaptation of the proof of Theorem 15.31 gives 


Theorem 5.4. Assume that each function Fi satisfies (|5.y^D -/E2)-/ [5. 6\) . and 
that H5.38\) - H5.39\) - /i5.40\ ) are satisfied for i = 1,... ,1, then 11^5.37^ admits at 
least one weak solution (pi, ..., pi) with each pi nonnegative. 
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6 Uniqueness 

We end the paper by a general uniqueness argument based on geodesic con¬ 
vexity. For the purpose of our paper it is enough to consider an internal 
energy that is a functional dehned on the subset Vs.cO^'^) of V{T'^) (consist¬ 
ing of absolutely continuous with respect to the Lebesgue measure on 
elements of V{T‘^)) and given by 

£{p) := [ E{p{x)) dx, pevu^^), 

JTd- 

where E : M+ —)■ M is a convex function with F'(O) =0 and E smooth on 
(0, -|-oo). Given p and u in Pac(T‘^), there is a unique optimal transport map 
T between p and p (see Cordero-Erausquin [5] and McCann [T3] for the case 
of a general Riemannian manifold) i.e. a map such that T^p = v and 

Wlip-,^) = f (f(T{x),x)p{x)dx. 

Moreover, this map is given by the gradient of a convex function T{x) = 
V0(x) with 0 convex on and such that x V4>{x) — x is periodic on 
(so that T indeed dehnes a map from T'^ to itself). In fact, on T^, we should 
rather write 

T{x) = exp 3 ,(—V'0(a;)) where 'i/j{x) = ~ 

Note in particular that d{x,T{x)) = |V'0(a:)| and then 

W 2 {p,u)= [ \V'ilj{x)\‘^p{x)dx. 

The Wasserstein geodesic between p and p is then the curve t G [0,1] i-G- Pt 
given by the McCann’s interpolation: 

vt := Tt#p, Tt{x) = exp,^(-fV0(a;)), 

it is indeed a constant speed geodesic: 

W2{yt, ^s) = \t- s\W2{p, y) 

and Tt is the optimal transport map from p to Vf Then the internal energy 
£ is said to be displacement convex whenever for every p and z/ in Pac(T‘’*), 
dehning Vt as above one has 

t G [0,1] I—)■ £{nt) is convex. 
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If E satisfies McCann’s condition: 


r E (0, +cxd) hE- r^E{r is convex nonincreasing 


( 6 , 1 ) 


then it is well-known that £ is displacement convex (see McCann ca) and 
that for p, -0 and v as above, one has (setting E\p) = pE\p) — E{p))\ 


£(u) — £(p) > lim 

i-)-0+ 


= lim 

t^o+ 



E{ut{x)) - E{p{x)) 
t 



VE'{p{x)) ■ Vip dx 



VE'{p{x)) ■ V'ipp{x) dx 


as soon VE'{p) in L‘^{p). Similarly, letting S {S{y) = y — V9{y)) be the 
optimal map from utop and using the fact that S{T{x)) = x i.e. V9{T{x)) = 
T{x) — X = —V'0(a)), we get 


£{p) -£{iy) > 


ijd 


VE\u{y))-V9{y)u{y) dy 


= - VE'{u{T{x)) ■ V9{T{x))p{x) da; 

= / '\/E'{iy{T{x)) ■ '\/'i/j{x)p{x) dx 


so that summing the two inequalities above gives 


0 > / (VE'(z/(T(x)) - VE'{p{x)) ■ \/ipix)p{x) dx. (6.2) 

Now let us consider the system for the evolution of I densities p = 

aip, = div(p,(Vi?'(A) + C,[p])), (6.3) 

on (0, -|-oo) X T'^. Let us assume that each Ei is a convex function as above, 
that Vi maps V{T^)^ into hL^’°°(T“') and that for some constant C one has 

|C,[p](a;) - V[p]{y)\ < Cd{x,y), 'ip E ViT^)^ and {x,y) eT^xT^ (6.4) 

and 

i 

iivcjp] - vv[u]\\L-iT^) v(p,z/) 6 nry x nry. 

i=i 

(6.5) 

Then the following uniqueness result holds: 
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Theorem 6.1. Assume that the drifts Vi’s satisfy ^6.4\ )- f67S\) and that the 
Ei’s satisfy McCann’s condition lid. 1\) . Let p and v be two solutions on 
(0,T) X of ad. A) such that 



with 

= 'V + Vi[ut\ 

then for every t G [0,T] one has the following stability estimate 


i 

^ ^ IT 2 {pi^tj 

i=l 


< e 


Act 


W^2^(Pi,0, ^i,o), 


2=1 


( 6 . 6 ) 

(6.7) 


( 6 . 8 ) 


which in particular gives unigueness for the Cauchy problem for \d. 31) . 

Proof. Using Theorem 5.24 and Corollary 5.25 from [20], assumption fl6.6p 
guarantees that 1 1 —)■ Wf{pi^t, V,t) is differentiable for a.e t G (0,T) with 


d 1 
dt 2 





Vi^t{.x))pi^t{,x)<lx 


where Vifi^t is such that Ti^t '■= exp 2 .(—VV'i.t) is the optimal map between p^ 
and Ui^f Now thanks to McCann’s condition and recalling fl6.2p . we have 


f Vifi^tix) ■ {VE-{ni^tiTi^t{x))) - VE'i{pi^t){x)))pi,t{x)dx < 0 , 

Jfd 

SO that, using 06.4^ - 06.51) yields 


<! 

[> 

■{Vi[ut]{At\ 

(x)) - Vi[pt\{x))pi,t{ 

x)da; 

JTd 






<1 

■{VHiAd 

A)) - Vi[ut]{x))pi^t( 

x)dx 

Jjd- 





+ f 


{Vi[vt]{x) - 

- Vi[pt]{x))pi^t{x)dx 


jTd 





< c 

[ WiA 

x)\d{x,Ti^t{- 

x))pi^t{x)dx 


J 

'jd 

r 


s 1 ' 


+ c[ 

/ 

{x)\‘^Pi^t{x)dx) ' V W2ipj,t, 

,t) 




i=i 



— C* I W^ipi^t, ^i,t) + Vd2{pi,t^ ^i,t) hU2(pyt, Uj^t) 

V j=i 
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Summing over i and using Gronwall’s lemma gives the desired inequality 

□ 

Remark 6.2. In the uniqueness result stated in Theorem 16.11 the integrability 
condition fl6.6p is made as an assumption since it ensures absolute continuity 
for the Wasserstein distance of the curves p and v. In the potential case where 
Vi = Vf/j, under the assumptions of Theorem 14.11 it can be checked (see [TT] 1 
that solutions constructed by the semi-implicit JKO scheme actually satisfy 

dSSl). 
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